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We review random loop representations for the spin-| quantum Heisenberg models, that 
are due to Toth (ferromagnet) and Aizenman— Nachtergaele (antiferromagnet). These 
representations can be extended to models that interpolate between the two Heisenberg 
models, such as the quantum XY model. We discuss the relations between long-range or- 
der of the quantum spins and the size of the loops. Finally, we describe conjectures about 
the joint distribution of the lengths of macroscopic loops, and of symmetry breaking. 
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1. Introduction 

Many electronic properties of condensed matter systems can be described by quan- 
tum Heisenberg models. Electrons are assumed to be localized but their spins in- 
teract with those of neighboring electrons. The study of quantum spin systems is 
notoriously difficult, and although several important results have been obtained, 
much remains to be understood. In this article we review probabilistic representa- 
tions. Probabilistic methods are not the most commonly employed in this subject, 
but they have proved useful nonetheless. A random walk representation has allowed 
Conlon and Solovej to obtain a lower bound on the free energy of the systemP This 
result was then improved by Toth with the help of the random loop representation 
described below Pa Another loop representation, this time for the Heisenberg antifer- 
romagnet, was proposed by Aizenman and NachtergaeleP It allows to relate the ID 
quantum model to 2D classical random cluster and Potts models. Recently, much 
progress has been achieved for another quantum model using similar probabilistic 
representations, the Ising model in transverse fields jHI l 16 l 20 l 

Random loop representations are also attractive per se. They allow to formulate 
open questions in a probabilistic setting, so that a whole new group of mathemati- 
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cians can reflect upon them. One should hope that probabilistic methods can shed 
a light on several properties of quantum spin systems. The situation now is rather 
the opposite: Several results have been obtained for Heisenberg models (absence 
of spontaneous magnetization in one and two dimensions, and its occurrence in di- 
mensions greater than two) using a genial combination of algebra and analysis, with 
insights from mathematical-physics. They provide stunning results when translated 
in the language of random loops. 

The quantum Heinseberg models and their probabilistic representations are in- 
troduced in Section [2j We describe in Section [3] the theorem of Mermin and Wagner 
about the absence of spontaneous magnetization/macroscopic loops in dimensions 
one and two, and the theorem of Dyson, Lieb, and Simon about the occurrence 
of spontaneous magnetization in greater dimensions. Sections [2] and [3] are mathe- 
matically rigorous. In the last Section |4] we describe the heuristics about the joint 
distribution of the lengths of macroscopic loops and we show that these conjectures 
are compatible with other conjectures concerning the breaking of U(l) or SO(3) 
symmetries. The claims of Section |4] are not proved. 

2. Quantum Heisenberg models and random loop representations 
2.1. Family of Heisenberg models 

Let (A, £ ) be a finite graph (simple, no loops), where A denotes the set of vertices, 
and £ denotes the set of edges. We consider the Hilbert space Ha = ® X £kH x , where 
each H x is a copy of C 2 . Let ti £ [-1,1] be a parameter. We consider the following 
family of Hamiltonians: 

' = —2 (^S^Sy + US^Sy + SpSy^j. ^ 

{x,y}e£ 

Here, S], = S^Id^^}, where S 1 , S 2 , S 3 denote the usual spin operators for spin-i 
systems (Pauli matrices). The partition function of the model at inverse temperature 
P is 

ZM(P,A) =Tre- /?ff A" ) . (2) 

We use the usual notation (•) to denote expectation of operators with respect to 
the Gibbs state, i.e. 

(A) = , } ln „ Trie-<'. (3) 

This family of Hamiltonians contains several cases of interest. 

• The case u = 1 gives the Heisenberg ferro magnet. In order to under- 
stand the physical motivation of this model, consider the Hilbert space 
for two spins at nearest- neighbors x and y, and the symmetry group of 
spin rotations. This group yields the irreducible decomposition T-l x <8> H y = 
singlet triplet. We want an interaction operator that is rotation invariant, 
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so it must be of the form CiP s j ng i e t + C2-Ptripict- Up to constants and a shift 
by the identity operator, we get ±S X ■ S y , hence the two Heisenberg models. 

• Choosing u = — 1, we get a model that is unitarily equivalent to the stan- 
dard Heisenberg antiferromagnet if the lattice is bipartite. The corre- 
sponding unitary operation consists of rotating all the spins of a sublattice 
by the angle 7r around the second spin direction. 

• The case u = gives the XY model. The standard representation involves 
interactions between spins in the axis 1 and 2, but the present choice is more 
suitable to the loop representation. It is well-known that the XY model is 
equivalent to the model of hard-core bosons, see e.g. Refs 111211 

2.2. Poisson point processes and loops 

The loop representation applies to general graphs (A, £). We consider a Poisson 
point process on £ x [0, /3] where 

• crossings occur with intensity "Mp. 

• bars occur with intensity ^-^ L - 

That is, the probability that a crossing occurs at edge e € £ and in the interval [£, t + 

e] is equal to ^^e + 0(e 2 ). Occurrences in disjoint intervals are independent events. 

(u) 

Let pg A denote the corresponding measure. Given a realization of this process, we 
define the loops in a natural way by following vertical lines, jumping to neighbors 
whenever a crossing or bars occurs (and continuing in the same direction in the 
case of a crossing, in the opposite direction in the case of bars). See Fig. [T]for an 
illustration. 




A 



Fig. 1. A realization of the two Poisson processes on edges and their associated loop configuration 
uonfx [0, /3] with periodic boundary conditions along the vertical direction. The number of loops, 
|£(u)|, is 6 here. This picture is for (A,£) being a line graph of 14 vertices, but arbitrary graphs 
can be considered. 
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Given a realization uj of the process p^\, let C(lo) denote the set of loops, and 
|£(w)| the number of loops. The relevant probability measure is 

^— 2^dp^\(u) (4) 

where the normalization is 

ZWC8,A) = J 2l £ Hl dp W(a,). (5) 

The case u = 1 and without the term 2> J ~'( U ^ is often called the random inter- 
change model. The size of its loops has been the object of several recent studies 
when the graph is a tre d 4 * 18 * 1 ^ and the complete graphPsEHU 



2.3. Relations between Heisenberg models and random loops 

A major advantage of the loop model is that it represents the Heisenberg models 
faithfully. Not only are the partition functions identical, but the spin correlations 
are given by the natural correlation functions in the loop model, namely that sites 
belong to the same loop. The following theorem has been progressively proved in 
RefsEEffi] 

Theorem 2.1. 

(a) The partition functions of quantum spins and random loops are identical: 

Z (M) (/3,A) = Z( u) (/3,A). 

(b) Spin correlations in directions 1 and 3 are given by 

{SlSl) = (SlSy) = ip^((x, 0) and {y, 0) belong to the same loop) . 

(c) Spin correlations in the direction 2 also have a loop counterpart: 

Here, E+ y denotes the event where (x, 0) and (y, 0) belong to the same loop, 
and the loop is moving in the same vertical direction at those points. E~ y is the 
similar event where the loop is moving in opposite vertical direction. Notice that 
E+ y UE~ y is the event where (x, 0) and (y, 0) belong to the same loop. An immediate 
consequence of the theorem is that 

\(S 2 x S%)\<(S 1 x S y :} = (S 3 x S 3 y ), (6) 

with equality if and only if u = ± 1 . 

Macroscopic loops are related to two physical properties of the system, namely 
spontaneous magnetization and magnetic susceptibility. This is stated in the follow- 
ing theorem. A more precise result, and a proof, can be found in Ref. 1301 Let £(o,o) 
denote the length of the loop that contains the point (0, 0) £ A x [0, 0\. 
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Theorem 2.2. 

(a) Relation between macroscopic loops and magnetic susceptibility: 

0' A V0|A|J |A| 2 g>77 2 g „=o 

(b) Relation between macroscopic loops and spin correlation functions: 

4 ST /Q3q3\ / 2d(l-tt) / ^(0,0) \ 4 ^ 3 3 

1 1 i,jeA V 1 1 11 x,y£A 

3. Results about spontaneous magnetization/macroscopic loops 

There are two major results, namely the absence of spontaneous magnetization in 
dimensions 1 and 2, and its occurrence in dimensions 3 and more, if the temperature 
is low enough. 

The first result is originally due to Mermin and Wagner and to Fisher and 
Jasnow for the decay of correlations.^ Their methods used Bogolubov's inequality 
and the translation invariancc in Z 1 or Z 2 . But the latter property is not essential. 
The theorem presented here and its proof are inspired by Frohlich and PfisterJ^ 
and NachtergaelePSlt is actually less general than those of Refs ll3|26l but its proof 
is simpler and it avoids the theory of abstract KMS states. 

Given a connected graph (A, £), we let d(x, y) denote the length of the minimal 
connected path from x to y. 

Theorem 3.1. Assume that there exists a constant C such that for any x £ A and 

any integer k, 

#{y:d(x,y) = k}<Ck. 

Then there exists a constant K , that depends on C and (3 but that does not depend 
on (A, £), such that for all i,i/£ A, 

< {slsl) < - = 

Vlogd(x,y) 

Notice that the bound also applies to correlations in other spin directions, using 
Eq. ©. It is not hard to check that the theorem rules out the possibility of spon- 
taneous magnetization on Z 2 and in other two-dimensional graphs, at arbitrary 
positive temperatures. 

Sketch proof. Let n = d(x,y). We consider the following unitary operation that 
rotates the spin at x but not at y, and that interpolates smoothly between them: 

U=Y[e^ s >, (7) 

zeA 
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with 

I otherwise. 
Since U*S%U = -Si and U*S 3 U = S 3 , we have 

Tr S 3 X S 3 e-P< } = -Tr S*S* e~^ H ^ u . (9) 

If we could remove the ?7's in the right side, the correlation would be zero. We show 
that it can indeed be done, at least approximately. We have 

U*H^U = - J2 U*(SlSl,+uS 2 z Sl+S 3 S 3 ,)U 

{z,z'}££ 

= ~ E e-»'& + ^(sle-*+*-+>W'Sl,^+*-+') !i ''' (10) 
{z,z'}es 

+ uSlSl + S 3 Z e-^'-W 2 ' S 3 Z , ) e -iMSl+sl). 

Observe that e~ laS i' S^, e laS i' = S^, — aS 3 , + 0(a 2 ), and a similar identity for the 
rotation of S 3 ,. Using the invariance of S^S^, + S 3 S 3 , and S^S 3 , - S 3 Sl, under 
rotations around the second direction, we get 

U*H^U = H^ + (^-cbz)(S 1 z S 3 ,-S 3 z Sl,)+ Od^-^f). (11) 

{z,z'}e£ {z,z'}££ 

We now use the identity Tr^4e~ H - Tr c~ H = lim s ^ ^Tr e ~ H + sA and Klein's 
inequality, so as to get 

- fl L—TrSlS 3 ( e-tzP - e^ u ' H ^ u ) = 
Z(")(/3,A) x y \ ) 

= ^4A) ^(e-<'^^ - e-^'W.a; ) + (s) (12) 

< sZ H (/3 ,A) Tr )U ~ O e"^^ + O(s). 

We use Eq. pT|) for the difference inside the trace. The middle term in Eq. (fTTj) 
gives because of the symmetry Y\ zeA c 17rS =, that sends onto — S*. We have 

i , ,i const ,„ . 

d(x, z) log(n + 1) 

for d(x, z) < n, otherwise. Using the assumption of the theorem, the higher order 
correction in (|11[) is bounded by 



E-L UU11SU . . 

k—— — - — -r < - — -. — . (14) 
(fclog(n + l)) 2 " Iog(n+l) 
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Combining Eqs ©, (fI2j). and (HH), we get for any s € (0, 1] 

- \ x »/ 2Z('")(/3,A) * "V / 



const , 
~ slog(n +1) 



(15) 



We choose s = , 1 in Eq. (IT21) and we get the claim of the theorem. n 

The second theorem is a positive result about the occurrence of spontaneous 
magnetization and it was proposed by Dyson, Lieb, and SimonfS] It was a "break- 
through based on another breakthrough" as Nachtergaele wrote itP^ It extended 
to quantum systems the method of infrared bounds developed by Frohlich, Simon, 
and Spencer for the classical Heisenbcrg modelP^l The result was initially proved 
for spin-i systems for d > 5, but it was extended to d > 3 by Kennedy, Lieb, and 
ShastrypS following observations by Neves and Perez Notice that it applies only 
to the antiferromagnetic case u € [—1,0]. It is a notoriously unsolved problem to 
extend this result to ferromagnetic systems. The theorem is elegantly formulated in 
terms of macroscopic loops. 

Theorem 3.2. Let (A, £) be a cubic box in Z d with d > 3, even side lengths, and 
periodic boundary conditions. Let u G [—1,0]. Then there exist (3q < oo and r\ > 0, 
independent of the size of the box, such that 



^< A V p\A\ 



E&0 (±m) > r,. 



It is possible to prove this theorem directly in the random loop setting without 
referring to the quantum frameworkP^l 



4. Heuristics and conjectures 

4.1. Joint distribution of macroscopic loops 

This section is based on the discussion of Ref. [15] Recall that a partition of the 
interval [0, 1] is a sequence (Ai, A2, A3, . . . ) of nonnegative, decreasing numbers such 
that J2i = 1- Let L\{uS) , ^(w) denote the lengths of the loops of uj in decreasing 
order. Then ( , j0t > ■ • ■ ) i s a random partition of [0, 1]. As |A| — > 00, we know 
that does not go to if d > 3 and f3 is large enough (see Thcorcm l3.2p . On the 

other hand, we should expect that a fraction of the domain A x [0, /3] belongs to loops 
of length /3, and more generally, to loops of bounded lengths. The first conjecture 
is a strong law of large numbers, that states that only finite and macroscopic loops 
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are present in the system. 

Conjecture 4.1. There exists v <G [0,1] such that 

lim lira ^. ) = 1 — v. ( finite loops ) 

lim lim > — tttt- = v. {macroscopic loops) 

k^oo |A|-kjo^ j3\A\ 1 y 

We know that f = in d = 1,2, and at high temperature in d > 3. But we 
should have v > if d > 3 and if the temperature is low enough. The conjecture 
also implies that ( %aT > ' ' ' conver g es hi distribution to a random partition 
of [0,1/]. 

In order to formulate a conjecture for the limiting distribution, we need to recall 
the definition of the Poisson-Dirichlet (PD) distribution. This is best done with 
the help of the closely related Griffiths-Engen-McCloskey (GEM) distribution. Let 
Xi,X2, ... be i.i.d. random variables of law beta(6*) (that is, Xi takes values in [0, 1] 
and F(Xi > s) = (1 — s) 6 for < s < 1). Then the random sequence 

X x , (1 - Xr)X 2 , (1 - X x )(l - X 2 )X 3 , . . . ) (16) 

is distributed according to GEM(6'). This is the "stick breaking" construction, since 
(1 — Xi) ... (1 — Xk) is what is left of the interval after chopping off k pieces. 
Rearranging these numbers in decreasing order, we get a random partition with 
distribution PD(6>). 

Conjecture 4.2. 

(a) If u = ±1, ( 

|A> ' L 0\A\l ' • • • ) converges in distribution to Poisson- 

Dmchlet(2). 

(b) If — 1 < u < 1, ( ' L /3\A\l ' ' ' converges in distribution to Poisson- 
Dirichlet(l). 

The mechanism behind the Poisson-Dirichlet distributions of the lengths of 
macroscopic loops is indirect but very general. The explanation is motivated by 
Schramm's work on the composition of random transpositions in the complete 
graphpS proving a conjecture of Aldous. This is explained in details in Ref. [15] 
in the case u = ±1, and it proceeds as follows. 

• Introduce a stochastic process such that the equilibrium measure 
2^ c ^^dpp l \(uj) is the invariant measure. 

• This yields an effective split-merge process on partitions. 

• The invariant measure of the split-merge process is Poisson-Dirichlet P^l 



The Markov process is quite natural and is defined as follows. 
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• A new edge-time (e, t) appears at rate y/2dt if its appearance causes a loop 
to split, and at rate l/\/2dt if it causes two loops to merge. 

• An edge-time already present disappears at rate if its removal causes a 
loop to split, and at rate l/y/2 if it causes two loops to merge. 

By considering all possible cases, we can check the detailed balance condi- 
tion: 

p(dw)2l £ (")lp(w, do/) = p(dw')2 |£( "' )l p(^', da;). (17) 

and since the process is ergodic, the measure 2^( u ^\ pp\(doj) is the unique invariant 
measure (up to a normalization). 

There is an important distinction between the case u = ±1 on the one hand, 
and the case —1 < u < 1 on the other hand. In the first case, any local change 
results in a merge or in a split (see Ref . fT5l for details). But in the case — 1 < u < 1, 
where both crossings and bars are present, it is possible that local changes do not 
break the loop, but only modify its internal order (think of o 8) The probability 
of splitting loops is then halved, and the effective split-merge process has stationary 
distribution PD(1) instead of PD(2). 

Models of spatial permutations are closely related. The occurrence of the 
Poisson-Dirichlct distribution can be proved in the "annealed" model where po- 
sitions are averaged uponP The lattice model is harder to study rigorously, but the 
mechanisms described above have been verified numerically.^ 

4.2. Macroscopic loops vs symmetry breaking 

Heisenberg models have natural rotation symmetries, namely SO(3) in the case 
u = ±1 and U(l) in the case — 1 < u < 1. They are expected to be broken in 
d > 3 and at temperatures low enough. These symmetries are not apparent in 
the loop representations. But we show here that there is a good reason why two 
different Poisson-Dirichlet distributions appear, PD(2) and PD(1), in these different 
situations. 

We consider the two-point correlation function (S x Sy) with x,y far apart. If 
u = +1, we expect that 

3(S 3 x S 3 y ) = (S x -S v ) = ^ JjS x ■ S y ) n dti, (18) 

where (-)^ is the pure state obtained by adding the external magnetic field efl ■ 
^2, x S x and by letting e \, after taking the thermodynamic limit. Using rotation 
invariance, it is enough to consider ft = e 3 . Then (S x Sy)g 3 = (S x Sy)g 3 = 0, and, 
as ||x — y\\ — >• oo, (S x Sy)g 3 = (S x > )g 3 (Sy)g 3 = with v the number appearing in 
Conjecture O Then (S%S$) = 

The case u = — 1 on bipartite lattices is similar. But one should perform the 
symmetry operation that gives the antiferromagnet + ^ S x ■ S y in order to use 
rotation invariance. We also get (S x Sy) = j^. 
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If — 1 < u < 1, we rather expect that 

2(SxSy) = (^x^y + Sx^y) = J ^x^y 



where il is in the plane (ei,e^). It is again enough to consider £1 
(S x S y }g 3 = (S 2 S 2 )g 3 = 0. Since (S%Sy)g 3 = (S^.)g 3 (Sy)g 3 = ^- as 



y\ 



(19) 

S3, and 
oo, we 



?3 o3 



find (S^, - T . 

Let us now calculate these correlations in the random loop representation. For 
x, y very far apart, the probability that they belong to the same loop is the same as 
the probability that two points in [0, 1] belong to the same element of the partition, 
multiplied by v 2 . This probability is easy to compute using the GEM distribution. 
Recall that if X is a beta(0) random variable, the expectation of X 2 is equal to 
(9+i)%+2 ) ' and tnc expectation of (1 — X) 2 is equal to -g^-. Using the definition 
(|16p for the GEM random partition, and summing over the probabilities that the 
numbers both belong to the kth element, we find 



[s, t belong to the same element) = V E((l - Xi) 2 ...(1-X k 



k>l 

E 

fe>i 

i 



9 + 2 



fe-i 



(0 + l)(0 + 2) 



(20) 



This gives 



0) and (y, 0) belong to the same loop) 




(21) 



This is indeed equal to 4(S x Sy), as found above. 

These heuristics calculations show that the conjectures about the Poisson- 
Dirichlet distributions of the lengths of the loops are compatible with the conjectures 
about the breaking of rotation invariance SO(3) or U(l). 



References 

1. M. Aizenman, E.H. Lieb, R. Seiringer, JR. Solovej, J. Yngvason, Bose- Einstein quan- 
tum phase transition in an optical lattice model, Phys. Rev. A, 70, 023612 (2004) 

2. M. Aizenman, B. Nachtergaele, Geometric aspects of quantum spin states, Comm. 
Math. Phys., 164, 17-63 (1994) 

3. G. Alon, G. Kozma, The probability of long cycles in the interchange processes, 
|http://arxiv.org/abs/1009.372"3"1 

4. O. Angel, Random infinite permutations and the cyclic time random walk, Discrete 
Math. Theor. Comput. Sci. Proc., 9-16 (2003) 

5. N. Berestycki, Emergence of giant cycles and slowdown transition in random transpo- 
sitions and k-cycles, Electr. J. Probab. 16, 152-173 (2011) 

6. N. Berestycki, G. Kozma Cycle structure of the interchange process and representation 
theory, |http://arxiv.org/abs/1205.4753| 



November 20, 2012 2:43 



WSPC - Proceedings Trim Size: 9.75in x 6.5in ueltschi-icmpl2 



11 



7. V. Betz, D. Ueltschi, Spatial random permutations and Poisson-Dirichlet law of cycle 
lengths, Electr. J. Probab. 16, 1173-1192 (2011) 

8. J. Bjornberg, Infrared bounds and mean-field behaviour in the quantum Ising model, 
|http://arxiv.org/abs/1205.3385| (2012) 

9. J. Conlon, J. P. Solovej, Upper bound on the free energy of the spin 1/2 Heisenberg 
ferromagnet, Lett. Math. Phys. 23, 223-231 (1991) 

10. N. Crawford, D. Ioffe, Random current representation for transverse field Ising model, 
Comm. Math. Phys. 296, 447-474 (2010) 

11. F.J. Dyson, E.H. Lieb, B. Simon, Phase transitions in quantum spin systems with 
isotropic and nonisotropic interactions, J. Statist. Phys. 18, 335-383 (1978) 

12. M.E. Fisher, D. Jasnow, Decay of order in isotropic systems of restricted dimension- 
ality. II. Spin systems, Phys. Rev. B 3, 907-924 (1971) 

13. J. Frohlich, C.-E. Pfister, On the absence of spontaneous symmetry breaking and of 
crystalline ordering in two-dimensional systems, Comm. Math. Phys. 81, 277-298 
(1981) 

14. J. Frohlich, B. Simon, T. Spencer, Infrared bounds, phase transitions and continuous 
symmetry breaking, Comm. Math. Phys. 50, 79-95 (1976) 

15. C. Goldschmidt, D. Ueltschi, P. Windridge, Quantum Heisenberg models and their 
probabilistic representations, in Entropy and the Quantum II, Contemp. Math. 552, 
177-224 (2011); [http://arxiv.org/abs/1104.0983l 

16. G.R. Grimmett, Space-time percolation, in In and out of equilibrium 2, Progr. Probab. 
60, 305-320 (2008) 

17. S. Grosskinsky, A. A. Lovisolo, D. Ueltschi, Lattice permutations and Poisson-Dirichlet 
distribution of cycle lengths, J. Statist. Phys. 146, 1105-1121 (2012) 

18. A. Hammond, Sharp phase transition in the random stirring model on trees, 



http://arxiv.org/abs/1202.1322l 



19. A. Hammond, Infinite cycles in the random stirring model on trees, 
http://arxiv.org/abs/1202.1319l 



20. D. Ioffe, Stochastic geometry of classical and quantum Ising models, in Methods of 
contemporary mathematical statistical physics, Lect. Notes Math. 1970, 87-127 (2009) 

21. T. Kennedy, E. H. Lieb, B. S. Shastry, Existence of N eel order in some spin-^ Heisen- 
berg antiferromagnets, J. Statist. Phys. 53, 1019-1030 (1988) 

22. T. Kennedy, E.H. Lieb, B.S. Shastry, The XY model has long-range order for all spins 
and all dimensions greater than one, Phys. Rev. Lett. 61, 2582-2584 (1988) 

23. E. Mayer- Wolf, O. Zeitouni, M.P.W. Zerner, Asymptotics of certain coagulation- frag- 
mentation processes and invariant Poisson-Dirichlet measures, Electron. J. Probab. 
7, 1-25 (2002) 

24. N.D. Mermin, H. Wagner, Absence of ferromagnetism or antiferromagnetism in one- 
or two-dimensional isotropic Heisenberg models, Phys. Rev. Lett. 17, 1133-1136 (1966) 

25. B. Nachtergaele, Quantum spin systems after DLS1978, in Spectral Theory and Math- 
ematical Physics: A Festschrift in Honor of Barry Simon's 60th Birthday, Proc. Sym- 



posia Pure Math. 76, 47-68, AMS (2007); http://arxiv.org/abs/math-ph/0603017 



26. B. Nachtergaele, unpublished lecture notes (2012) 

27. E.J. Neves, J.F. Perez, Long range order in the ground state of two-dimensional anti- 
ferromagnets, Phys. Lett. A 114, 331-333 (1986) 

28. O. Schramm, Compositions of random transpositions, Israel J. Math. 147, 221-243 
(2005) 

29. B. Toth, Improved lower bound on the thermodynamic pressure of the spin 1/2 Heisen- 
berg ferromagnet, Lett. Math. Phys. 28, 75-84 (1993) 

30. D. Ueltschi, Random loop representations for quantum spin systems, in preparation 



